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ABSTRACT 

We study the spherical collapse model for several dark energy scenarios using the fully nonlin- 
ear differential equation for the evolution of the density contrast within homogeneous spher- 
ical overdensities derived from Newtonian hydrodynamics. While mathematically equivalent 
to the more common approach based on the differential equation for the radius of the pertur- 
bation, this approach has substantial conceptual as well as numerical advantages. Among the 
most important are that no singularities at early times appear, which avoids numerical prob- 
lems in particular in applications to cosmologies with dynamical and early dark energy, and 
that the assumption of time-reversal symmetry can easily be dropped where it is not strictly 
satisfied. We use this approach to derive the two parameters characterising the spherical- 
collapse model, i.e. the linear density threshold for collapse 6^ and the virial overdensity Ay, 
for a broad variety of dark-energy models and to reconsider these parameters in cosmologies 
with early dark energy. We find that, independently of the model under investigation, ^c and 
Ay are always very close to the values obtained for the standard ACDM model, arguing that 
the abundance of and the mean density within non-linear structures are quite insensitive to the 
differences between dark-energy cosmologies. Regarding early dark energy, we thus arrive at 
a different conclusion than some earlier papers, including one from our group, and we explain 
why. 
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1 INTRODUCTION 



• Over the last decade, a wealth of evidence was accumulated in 
favour of the conclusion that the expansion of our Universe is ac- 
, celerating , mainly fr om the obs ervation of the type-la supemovae 
■ JRiesset al. 1998; Perlmutteret al. 1999; Kowalski et al. 2008) and 
the cosmic microwave background (CMB) (Komatsu et al. 2009) in 
combination with mea surements of the Hubble constant and large- 
scale structures (LSS) JCole et al 120051) . Assuming the validity of 
general relativity on large scales, a possible explanation for this ac- 
celerating expansion, is obtained by introducing a component of 
the cosmic fluid, the dark energy, with equation-of-state parameter 
w < -1/3. 

Despite efforts from both observational and theoretical sides, 
the nature of the dark energy remains obscure. Consequently, 
a plethora of different models has been proposed for the ori- 
gin and the time evo lution of the dark energy, see for example 
ICopeland et al.l ( l200q) for a comprehensive review. The simplest 
model assumes that the dark energy is connected with the vacuum 
energy, the so called cosmological constant, with equation-of-state 
parameter w = —1. Despite the fact that observations constrain the 
present value of w quite tightly, the time evolution of the equation- 
of-state parameter is rather poorly constrained. Thus, it is natural to 



E-mail: francesco@ita.uni-heidelberg.de 



study more general classes of models allowing a time evolution of 
the dark-energy component, such as models involving scalar fields. 

Scalar fields occur naturally in particle physics and in string 
theory and could thus be candidates to explain the nature of the 
dark-energy if they are sufficiently strongly self-interacting. This 
class includes quintessence models, phantom models, K-essence, 
tachyon models and so forth. Scalar fields are described by their 
Lagrangian density with a kinetic term <^^/2 and a potential 
term V{cj>). The equation-of-state parameter then follows from the 
canonical energy-momentum tensor. If the dark energy is spatially 
homogeneous, P = (1)^/2 + V{(j)) and p — </>^/2 — V{4>), giving 
w = P/{pc^). 

The cosmological-constant case is recovered if the kinetic en- 
ergy is negligibly small compared to the potential energy. Dark en- 
ergy affects first of all the expansion rate, causing geometrical ef- 
fects that can be revealed through distance measurements, such as 
the luminosity distance to distant supemovae. Secondly, it affects 
structure formation, the early stages of which can be quantified by 
the growth factor. Thus, structure formation will be affected by the 
amount of dark energy and by its dynamical evolution over cosmic 
history. 

One model recently suggested as a candidate for solving 
the fine-tuning problem of the cosmological c onstant was the 
class of early dark-ener gy cosmologies (EDE) l lWetterichI 1 2004 ; 
iDoran & Robbers! 1 20061) . according to which the contribution of 
dark energy at early times is not negligible. Thus, to produce the 
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same amount of structure now, structure formation should start ear- 
lier and proceed more slowly than in the common ACDM model. 
It would thus compensate the additional, opposing effects of the 
early dark-energy contribution. An analytic calculation based on 
the evolution equation for the radius of a spherical, homogen eous 
perturbation and various assumptions JBartelmann et alj|2006n im- 
plied a substantial increase in the number of objects compared 
to a stan dard ACDM model. While this expectation was con- 
firmed bvlSadeh et al.l ( l2007h. subsequent A^ - body simulations by 
iFrancis et alj <2009ah and lOrossi & Springell 120091) found instead 
that the effect of EDE on the mass function of dark-matter haloes 
and its evolution is almost negligible: The EDE class of models 
predicts differences in the mass function of only a few percent with 
respect to the ACDM model. A new numerical derivation of the 
linear overdensity parameter, also ba sed on the differentia l equa- 
tions of the spherical collapse model JFrancis et alj|2009a) was in 
perfect agreement with the numerical simulations. 

Motivated in part by this discrepancy of results derived from 
the same model, we are here addressing the problem of determin- 
ing the time evolution of the linear overdensity 5c[z) in a com- 
pletely different way, using a perturbative approach based on New- 
tonian hydrodynamics directly. One of the advantages, apart from 
increased numerical stability, is that no time-reversal symmetry 
needs to be assumed for the spherical collapse. 

The structure of the paper is as follows. In Sect.[2l we present 
the basic equations of Newtonian hydrodynamics and sketch the 
derivation of the equations used to obtain the linear over density 
threshold Sc and the non-linear virial overdensity parameter Ay. In 
Sect. [3] we briefly describe and motivate the cosmological models 
investigated in this work, while we compare them in Sect. |4] with 
the ACDM model. We present our conclusions in Sect. (6) In the 
appendixlA] we discuss why the prev ious theoretical es t imatio ns of 
5c for the EDE m odels obtained by iBartelmann et al.l J2006l) and 
ISadeh et al . f2007^ differ from the results obtained in this work and 
bv lFrancisetal.t2009h) . 



2 NEWTONIAN HYDRODYNAMICS OF A 
RELATIVISTIC FLUID 

We review here the derivation of the differential equation determin- 
ing the evolution of an overdensity 5. The final non-linear equation 
specialised to w = is not new, but has already been used by 
several authors in the cont ext of structure formation ( Padmanabhan 
ll996l : lAbramo et alj2007h and for the study of the spherical a nd el- 
IJDSoidal coUapse teemard"eaulll994l0htaet al.ll2003ll2004h . The 
linearised equation was presented in lColes & LucchinI J2002i) spe- 
cialised for two limiting cas es, namely dust (w = 0) and relativis- 
tic matter {w = 1/3), and in lLima et alj ( ll997h for a general model 
with constant w. 

Our study, based on the work bv lAbramo et al.l <2007D where 
the equation for the evolution of the overdensity 5 was generalised 
to allow for time-dependent equation-of-state for the dark energy 
component, has two novel aspects: First, we generalise the evolu- 
tion equation to a generic collapse geometry, rendering the spher- 
ical and ellipsoidal models special cases. Second, the newly ob- 
tained generality of the method allows its application to modified- 
gravity cosmologies and coupled-quintessence models. This will be 
postponed to future work. 

Following the work by lAbramo et al.l 120071), we derive our 
equation including the pressure terms explicitly and assuming that 
the fluid satisfies the equation-of-state P = wp(? . Starting from 



this point, the final equation will be in its most general form and 
can then be specified to a particular fluid simply by adopting the 
appropriate equation-of-state. 

We start from the continuity equation for the energy- 
momentum tensor in general relativity, V,jT''^" = 0. For a perfect 
fluid, we have 



T"" = (pc^ + P)u''u'' + Pg^' 



(1) 



where p is the density of the fluid, P its pressure, u its 4- velocity 
and g'^'^ the metric. 

Contracting the continuity equation once with u^ and once 
with the projection operator g^a +u^Ua one obtains the relativistic 
expressions for the continuity and the Euler equations, respectively: 



ot c^ 

dv ,^ „ w ^ ^ c^VfP + vP „ 
ot pc' + P 



(2) 
(3) 



Here v is the velocity in three-space, $ is the Newtonian gravita- 
tional potential and r is the physical coordinate. 

The 0-0 component of Einstein's field equations gives the rel- 
ativistic Poisson equation 



V"$ = 47rG(p+-^ 



(4) 



The continuity equation for the mean background density, obtained 
from the spatial components of Einstein's equations, is now modi- 
fied to 



p+3H{p+ 







(5) 



where p = — ""'^ is the background mass density of all contri- 
butions to the cosmic fluid, and fifluid is its density parameter. 

As usual, we introduce comoving coordinates x = f/a and 
define 

(6) 
(V) 
(8) 
(9) 

where H(a) is the Hubble function and u{x, t) is the comoving 
peculiar velocity. Inserting Eqs. l6}-([9} into Eqs. ([^J-Q, we find 
the equations 



p{x, t) = 


= p{l + 5{x,t)), 


P(x,t) -- 


= wp{x,t)c^ , 


$(f,t) -- 


= $o{x,t) + (l>{x,t) , 


v{x, t) = 


- a\R(a)x ^ u(x,t)\ 



<5 + (l + w)(l + (5)Vs--u = 

^ + 2//^+(u- Vs)72+-t-V3(?!> = 
at a-' 

V|0 - 47rG(l + 'iw)a'p5 = 



(10) 
(11) 
(12) 



We now take the divergence of the Euler equation jilt and 
recall the decomposition 



V -[{u- V)u] = -0' + (J- - uj- 



(13) 



into the expansion — Va • u, the shear tensor a^ = aijcr^-' and 
the rotation tensor oj^ = ujijuj''-' . Next, taking the time derivative 
of the continuity equation UOb and combining all three equations. 
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we arrive at the fully non-linear evolution equation 
S+(2H ^ ^^ 4 + 3"' '" 



l + wj" 3{l+w)l + 5 
4nGp{l + w){l + 3w)5{l + <5)- 

{l + w){l + S){a^-LO^)^0 



(14) 



Note that the shear is a symmetric traceless tensor, while the rota- 
tion is antisymmetric. They read 



1 fdv^ &tl_ 



\05., , 



(15) 



Specialising Eg. ( [141 for dust {w = 0), we recover Eq. (41) of 
Ohta et al. ('2003'). We notice that Eq. (O generalises Eq. (7) of 
Abramo et al. (2007) to the case of a non-spherical configuration 
of a rotating fluid. 

Changing the independent variable from the time t to the scale 
factor a using the relation dt = aH{a)da, the evolution equation 
assumes the form 



a bj 1 



5'- 



4 + 3u> (5'2 
3(1 + io) 1 + 5' 



3 Jlfluid£_ 
2a2£2(a) 



g{a){l + w){l + ■iw)5{l + 5)- 

-(l + u.)(l + <5)((7'-a;") = 



(17) 



where Hfluid.o is the density parameter of the fluid at a = 1, and 
g{a) is a function specifying the time evolution of the dark-energy 
model considered. 

In the following, since we are interested in the collapse of a 
homogeneous sphere, we ignore the rotation and the shear tensors. 
The shear tensor vanishes for a sphere anyway. We will also restrict 
the treatment to spherical perturbation filled with dust, having w = 
and g{a) = a^^ . Thus, the non-linear equation to be solved reads 

3 f^m.O 



5(1 + 5) = . (18) 



\a E ) 31 + 5 2a5£;2(a) 

We notice that Eq. ( |17t has a singularity when w — —1. To see 
what happens for the cosmological constant, we multiply both sides 
with 1+w and then specialise tow = —1. We obtain 5'^/(l + 5) — 
0, implying 5 = const, and with appropriate initial conditions, the 
constant can be set to zero, so the cosmological constant can not 
clump as expected. 

2.1 Determination of 5c and Av 

The linearised Eq. dlSt reads 

and its solution, for appropriate initial conditions, will give the lin- 
ear overdensity parameter 5 at any point in time. This equation is 
also used to determine the time evolution of the growth factor if 
suitable initial conditions are used. 

To determine the appropriate initial conditions, we start by 
considering Eq. dlSl l. We know that, since this represents the non- 
linear evolution of the density contrast, its value at a some cho- 
sen collapse time diverges, 5 — >■ oo. Thus, we search for an initial 
density contrast such that the 5 solving the non-linear equation di- 
verges at the chosen collapse time. Numerically, we assume this to 



10^ 




8 


1 1 1 1 1 

5 linear 

non-linear 


'' : CD ' '' 

E 




-. 


103 


I- 






egim 
ear r 




-: 


10^ 


- 










- 


in^ 


r 






:i : g 




- 








5^=1.686 


< »- 




- 


10" 


- 








,--"'*'>•■ 












,-- CO ■ 




in-^ 


- 






^^'^ 


O ■ 


- 
















10-2 


I- 




_,^--' 




CO ; 


- 


10-^ 


r 


_,- 


....-'" 




8i 


- 


Hn-4 


■'' 




, , 1 , 


,1 : , , 1 


, , 1 


- 



10' 



10" 



lO" 



10" 



10" 



Figure 1. Linear (blue dashed curve) and non-linear (cyan short-dashed 
curve) evolution of the overdensity pai'ameter S. An EdS model and a sphere 
collapsing at z = are assumed. We notice how, after the initial parallel 
evolution, the non-linear solution grows very fast with the scale factor, in 
comparison to the linear solution. 



be achieved once 5 ^ 10 . Since the curve S{a) representing the 
non-linear density evolution turns very steep towards the collapse, 
the result is very insensitive to the exact choice of this threshold 
value as long as it is a large number. Once the initial overdensity is 
found, we use this value as an initial condition in Eq. ( |19l l to find 
5c. 

Since we are dealing with second-order equations, two initial 
values have to be given, one for the initial overdensity Si and the 
other for the initial rate of evolution, 5[. We know that initially, 5- 
for the sphere should be small, thus we set it to 5; = 5 x 10~^, 
corresponding to the initial scale factor used for starting the in- 
tegration of the two differential equations. We carried out several 
numerical tests to check the dependence of the solution on 5; and 
found that the result does not depend on the precise value of 5{. 
Specifically, we considered several values in the interval between 
10"'' and 10"* and saw perfect convergence of the solution. Also 
setting 5[ = 0, the result does not change considerably. 

In Fig.[T] we show the solution 5 of the non-linear (cyan short- 
dashed curve) and the linear (blue dashed curve) evolution equa- 
tions as a function of the scale factor, for an EdS model, supposing 
that the sphere collapses at z — 0. We see that the linear solution 
grows linearly with time, reaching the correct value for 5c = 1.686 
at the collapse scale factor, while, after developing in parallel ini- 
tially, the non-linear solution starts deviating and growing expo- 
nentially. To obtain the virial overdensity, supposing that dark en- 
ergy does not clump, we follow the pr escription o f Maor & Lahavl 
( l2005b , which generalises the work bv lWang &"steinhardll ( 11993) 7 

Knowing the non-linear time evolution, it is possible to infer 
all the other properties of the collapsing sphere, in particular the 
time evolution of the radius, the turn-around scale factor ata. when 
the sphere reaches its maximum radius, and the overdensity at turn- 
around (". The virial overdensity is defined as Av = S^i + 1 = 
C,{x/y)'^, where x = a/ata is the normalised scale factor and y is 
the radius of the sphere normalised to its value at the turn-around. 

To determine the turn-around scale factor, we solve Eq. ( |17| l 
and determine the quantity log(5ni + l)/a^. Apart from a multi- 
plicative constant, this is the inverse of the collapsing sphere's ra- 
dius and assumes a minimum at the turn-around scale factor ata. 
To determine the virial overdensity at turn-around C„ we integrate 
Eq. dni l up to ata and add the result to unity. 
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3 THE MODELS 

As mentioned before, the nature of the dark energy is still unknown, 
which leaves room for a plethora of phenomenological or ad hoc 
models being discussed in the literature. Here, we briefly review 
some models, characterised by the requirements that they try to 
explain the accelerated expansion of the Universe in terms of a 
smooth and slowly varying component, the dark energy, quanti- 
fied by a certain equation-of-state parameter w{a) and formulated 
in the framework of general relativity. The dark-energy component 
satisfies the following evolution equation 



Table 1. Parameter values for the quintessence models. 



P + 3//P+— =p + 3H(1 + w)p = 



(20) 



Our co llec tion of models i s based on the works bv lSzvdlowski et al.l 
( 120061) and ljennings et al.i ([2009J where A'^-body simulations of dif- 
ferent quintessence models are studied. We shall use the following 
cosmological parameters: Q.m,Q ~ 0.274, Jlq.o = 0.726. For flat 
models Hk^ = 0, while for models with a curvature term we set 
^K,o = -0.018. 

3.1 ACDM model 

The simplest model used to fit the data that explains the late-time 
accelerated expansion of the Universe has a cosmological constant, 
with equation-of-state parameter lu = — 1 independent of time. Be- 
cause of this, the contribution of the cosmological constant starts 
dominating only recently and becomes rapidly negligible towards 
higher redshift, such that, at high redshift, it converges towards the 
EdS model. Despite conceptual problems associated with it, it is 
currently the simplest model fitting virtually all available observa- 
tional data. We shall thus assume a spatially flat ACDM model as a 
reference model, since observations suggest negligible spatial cur- 
vature. However, we will also consider finite curvature, even if this 
parameter, according to the current limits, is quite small. 



Model 


wo 


Uim 


dm 


Am 


INVl 


-0.4 


-0.27 


0.18 


0.5 


INV2 


-0.79 


-0.67 


0.29 


0.4 


2EXP 


-1.0 


0.01 


0.19 


0.043 


AS 


-0.96 


-0.01 


0.53 


0.13 


CNR 


-1.0 


0.1 


0.15 


0.016 


SUGRA 


-0.82 


-0.18 


0.1 


0.7 



rolling down its potential reaching an attractor solution. An inter- 
esting feature is that the scalar field tracks the dominant component 
of the cosmic fluid. The second class instead keeps the ratio be- 
tween the density of the scalar field and that of the background con- 
stant. ThejnodelsJNVl_Mid INV2 h ave inverse power - law poten- 
tials JCorasaniti & Copeland 2003; S anchez et alj2009l : ICorasanit 



l2004f) . the SUGRA model is a t ypical example for t racking be- 
haviour as well as the CNR model I Copeland et alj2000,). while the 
mode ls 2EXP JBarreiro et ^l2000t) and AS JAlbrecht & Skordij 
1 200011 are examples for scaling fields. 

Given the appropriate equation-of-state parameter for each 
model, we can solve the continuity equation |20] leading to the so- 
lution 



Poe 



-3 J-^[l-\-iu{a )]dlna 



(21) 



Thus, the expansion function reads 



E{a) = 



^K,0 , o 

— 5 — hi2Q,oexp 



"i±ifM^., 



da> 



(22) 

All these models can be described by the following equation-of- 
state parameter, valid after matter-radiation equality. 



3.2 Quintessence models 

An immediate extension of the cosmological-constant scenario is 
described by a scalar field very weakly or not interacting with the 
matter component. This scalar field can be in principle the inflaton 
itself, even if the vast majority of the scenarios assumes it to be 
independent of the scalar field actually driving the observed accel- 
erated expansion. These models are described by a kinetic energy 
and a potential energy characterised by a given functional form, 
that can either be motivated by theory or introduced ad hoc, such 
as power-law potentials. The equivalent mass of the scalar field is 
given by the second derivative of the potential term. Compared to 
the cosmological-constant case, these models have a time-evolving 
equation-of-state parameter. They are justified by the fact that a 
time-dependent equation-of-state parameter, not excluded by ob- 
servations, naturally arises in the framework of a scalar field the- 
ory. It is assumed that dark energy does not clump, at least not on 
the relevant scales accessible to cosmological studies. The sound 
speed of quintessence models is directly related to the equation-of- 
state parameter by the relation Cs = ^/wc. 

Here, together with the early dark-energy models (see 
Sect. |3.2.Tt , we study the six models used by Jentiitigs et a],: (2 009i) 
to which we refer the reader for more detail. All models can be di- 
vided into two broad classes, tracking scalar fi e lds JSteinhardt et al.l 
19991) and sca ling fields faalliwelll 1 1981 IWands et alj 1 19931 : 
Wetterichlll995r) . Tracker fields are characterised by a scalar field 



w(a) = Wo + {wui — wo)- 



l + e'^ 



(23) 



1 + e" 



In Tab. 13.21 we list the values of the parameters am, Am, Wm 
and Too characterisin g the models discussed above. 

The CPL model IChevalUer & Polarskil ( l200lh : lLiride3 J2003h 
has the equation-of-state parameter 



7(z) = -l 



1 + Z 



(24) 



mimicking a cosmological constant at low redshift and growing up 
to zero at very early times. 

In Fig.|2]we show the equation-of-state parameter of as a func- 
tion of the scale factor a for the different models. The linestyle 
coding is given in the figure. We notice that these models show 
different behaviour: the INVl and INV2 models show a gentle in- 
crease of the equation-of-state parameter, which is almost constant, 
except for late times. The other models instead show a very large 
change in to at late times, reaching a constant value quite soon in 
cosmic history. Also, the values of the equation-of-state parameter 
at a = 1 cover a broad range, from to — —0.4 for the INVl model 
till TO = — 1 for essentially all others. 

Even if th e SUGRA model was s tudie d both analyti c ally an d 
numerically by iMainini et al] ( l2003bl) and iMainini et all ( l2003ah . 
we determine once again the expected Sc for this model as a test of 
the validity of our approach and for completeness. 



Spherical collapse model in dark energy cosmologies 5 




-1.4 



INV1 

INV2 

2EXP 

AS 

CNR 

CPL 

SUGRA 



10" 



10" 



10"' 



10"- 



Figure 2. Time evolution of the equation-of-state parameter as a function 
of tlie scale factor a for the different quintessence models studied in this 
work. The light-green dashed and the dark-green short-dashed curves rep- 
resent the INVl and INV2 models, respectively; the blue dotted curve the 
2EXP model, the cyan dot-dashed curve the AS model, the orange dot- 
short-dashed and red dot-dot-dashed curves the CNR and the CPL models, 
while the solid violet curve represents the SUGRA model. 



3.2.1 Early dark energy models 



A part icular class of dark-ener gy models, introduced bylWetterichl 
( I2OO4I) and studied in detail bv lDoran & RobbersI ( l200dh . introduc- 



ing a convenient functional form for its equation-of-state parame- 
ter, has a small but finite density of early dark energy (EDE). This 
class of model, where the density parameter of the dark-energy 
component r emains at the level of a fe w percent at very early times, 
was used by Bartelmaim et alj ( |2006[) for their study of non-linear 
structure formation. 

Its equation-of-state parameter is implicitly given by 



3w(a) — 



a + ac 



nci{a)il~nci{a)) = -^^^, (25) 



where Oeq is the scale factor at matter-radiation equality, and 
Qq (a) represents the time evolution of the dark energy component, 



fiQ(a) 



^Q,0 — f^Q,c(l 



-:iwo\'r 



f^Q,0 + f^n 



+ t2Q,e(l-a-^"'°)^ 



(26) 
where f2Q,o is the density parameter of dark energy today, f2Q,o its 
density parameter at early times, wo the present equation-of-state 
parameter, and 7 is a shape parameter controlling the importance of 
the terms containing J1q,c- The expression explicitly assumes a flat 
universe. We adopt 7 = 1 here. Since the equation-of-state param- 
eter for the EDE model is very similar to that of the quintessence 
models mentioned above, we do not report it here. Further detail 
on the comparison of the new approach with the old one to deter- 
mine 5c for the EDE models will be given in Appendix [A] where 
w e also compare our predic tion with the numerical mass function 
bv lGrossi & SpringellOOOgI) . 

Here, we use the cosmoloeical pa rameters of the model EDE4 
from lWaizmann & BartelmannI ( l2009h . i.e. fim.o = 0.282, JIq.o — 
0.718 and wo ~ —0.934, while for the comparison with the nu- 
merical simulations, we use f2m,o ~ 0.25, nQ,o ~ 0.75 and 
Wo = -0.99. 



3.3 Chaplygin gas and Casimir effect 

An alternative to scalar fields for explaining the current accel- 
eration proceeds by specifying an exotic equation of state sat- 
isfying the condition for acceleration, w < —1/3. One exam- 
ple is the Chaplygin gas, first proposed in aerodynamics and 
subsequently derived from the action in string theory JOgawal 
UOOCt) . It was use d in cosmology as a possib l e alternative t o dark - 
energy models by Kamenshchik et al.l <200lh : lFabris et alj ( l2002h : 
ISzvdlowski & Czaial J2004h . The equation of state of the gener- 
alised Chaplygin gas assumes the form 

— ^^ 

where A > and a are constants. The classical Chaplygin gas is 
recovered for a — 1. Using the continuity equation for the gener- 
alised Chaplygin gas, one obtains the dependence of the density on 
the scale factor 

i/(i+a) 



(27) 



P = 



A + 



B 



,3(0 + 1) 



(28) 



where B is an integration constant. The equation-of-state parameter 
can be written in the form 

A 



w{a) 



A + 



(29) 



„3(a + l) 

where A = -t(;o(nQ,o/9c)'+" and B = (1 + wo)(nQ,oPc)'+". 
/9c is the present critical density and wq — —A/{A + B) is the 
present value of the equation-of-state parameter. We study the clas- 
sical Chaplygin gas with a = 1 and a generalised version with 
Q — 0.2. Both models have wo = —0.8. 

Another possible way to explain the current accelerated ex- 
pansion is to study quantum properties of the vacuum using the 
Casimir effect. This effect arises from a change in the zero-point 
oscillation spectrum of a quantised field when the quantisation do- 
main is finite or the space topology is non-trivial. In a cosmological 
context, the Casimir effect is relevant if the topology is not simply 
connected or when compact extra dimensions are involved. In a 
more general setting, it can be used to study the properties of the 
vacuum. In this context the contribution given by the Casimir force 
is scaling like a^*, exactly like relativistic species do. 

The expansion function is 



E{a) 



fin 



+ ^Q,0 - 



ftc 



(30) 



where f^cass.o is the density of the Casimir component today. If 
one wants to interpret this as a time evolution of the dark energy 
component, one can invert Eq. OOt using the general equation 



w[a) — 



^_^2^dlnE(al 



3 a'-'Eia)'^ 



1 



n„ 



a^B(a)^ 



(31) 



(32) 



aaB(a)^ 

to obtain the equation-of-state parameter 

, , 1 Srio.oa + nCaas,0 

"^("^ " ~^Ti ^4 — a 

6 liQfia* — SiCass.O 

Here, we shall assume ficass.o ~ —0.00035. 

In Fig. [3] we show the time evolution of the Chaplygin (ma- 
genta dashed curve) and of the generalised Chaplygin (turquoise 
dotted curve) gas, and in brown dot-dashed the equation-of-state 
parameter for the Casimir effect. The two curves representing the 
generalised Chaplygin gas are very similar, only the initial slope 
changes with the change of the a parameter. 
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Figure 3. Time evolution of the equation-of-state parameter. Tiie (turquoise 
dotted) magenta daslied curve shows the (generalised) Chaplygin gas while 
the brown dot-dashed curve represents the model based on the Casimir ef- 
fect 



3.4 Phantom models and topological defects 

A better fit to type-la supernova data is achieved if one assumes 
a varying equation-of-state parameter with some phantom crossing 
at low redshift, or in other words that there w < ~1. Models ful- 
filling this condition are called phantom models, and they seriously 
challenge the foundations of theoretical physics, since they violate 
several energy conditions. Here, we study five different phantom 
models, all having a constant equation-of-state parameter, in partic- 
ular we focus on the models with w = —4/3, w = —3/2, w = —2, 
w — —3. Certain grand unified theories predict topological defects 
to have formed in the early universe and since they rapidly diluted, 
their abundance should be very low. Here, we consider models with 
w = -2/3. 



4 RESULTS FOR Sc AND Av 

In this section, we discuss the results for the linear overdensity 
parameter and the virial overdensity for the models introduced in 
Sect. [3] keeping the ACDM model as a reference. 

Our main results are shown in Fig.|4] The right panels show re- 
sults for the virial overdensity Av (z), while the left panels are spe- 
cialised to the linear overdensity 5c ■ The upper panels refer to the 
quintessence models, the middle panels refer to the (generalised) 
Chaplygin gas and to a cosmology with Casimir effect taken into 
account. The lower panels show results for the models with a con- 
stant equation-of-state parameter (non-flat ACDM model, topolog- 
ical defects and phantom models). Linestyle labels are explained in 
the figure caption. 

From Fig. [4] it is quite evident that all models considered, in- 
cluding the EDE cosmologies discussed in the appendix, behave 
very similarly to the flat ACDM cosmology, irrespective of the 
equation-of-state parameter, be it constant or varying with time. 
At 2; = 0, the difference in 5c is at most of 2% for the gener- 
alised Chaplygin gas, while for the very large majority it is even 
less. All models asymptotically approach the EdS limit at high red- 
shift. This result, that all models give essentially the same results, is 
quite important: It shows that the linear density threshold Sc from 
the ACDM model is very close to the precise value in other cos- 
mologies even if the equation-of-state parameter considerably dif- 
fers from w — —1. We argue that a possible enhancement in struc- 



ture formation might be caused by rapidly varying or discontinuous 
equation-of-state parameters, for example if they contain bumps or 
peaks. From a physical point of view, huge differences from the 
ACDM models might result from modified-gravity scenarios, such 
as coupled dark-energy models. 

It is also interesting to see that the equation-of-state parameter 
has very little impact on the evolution of 5c ■ We argue that this can 
be due to the fact that the equation-of-state parameter is always in- 
tegrated over and thus its effects are smoothed over cosmic history. 
It would be interesting to work out with the equations governing 
the evolution of the overdensity which conditions must be satisfied 
by the equation-of-state parameter to have significant effects on 5c. 

The same considerations apply to the virial overdensity Av. 
Deviations at low redshift are at most of the order of a few per- 
cent, thus having negligible impact on non-linear structure evolu- 
tion. This fact has also a practical advantage: all quantities depend- 
ing on Av will be virtually unaffected if the virial overdensity of 
the ACDM model is used as an approximation. 

We thus conclude that the models discussed above have no 
significant impact on non-linear structure formation, including 
the EDE models presented in the appendix. Hence, the conclu- 
sions on the thermal S uny aev-Ze l'dovich effect (Sadeh et al. 20()7|; 
Waizmann & Bartelm anrJ 120091) , lensing (IFed eli & BartelmannI 
,2007. : .Fedeli et al., ,2008h and clustering JFedeh et al., ,2009.) are 
based on erroneous assumptions. 



5 VOLUME EFFECTS ON HALOES NUMBER COUNTS 

In the previous section, we studied the impact of different dark 
energy models on the linear overdensity threshold 5c. A quantity 
closely related to observations is the mass function, representing 
the number of collapsed objects per unit mass and volume. Since 
it only depends on 5c and on the growth factor, no appreciable dif- 
ferences are expected between the models studied. An important 
quantity that can be derived from observations is the total num- 
ber of haloes above a given mass in a complete survey volume. 
The minimum mass detectable in a survey is generally a function 
of redshift and changes with the observed wave band; moreover, it 
will also depend on the survey according to the instrument sensitiv- 
ity. Since we do not intend to specify an individual survey here, we 
assume the minimum mass to be independent of redshift. An ideali- 
sation in this approach is that the catalogue of objects is considered 
to be complete in order to compare observations with theoretical 
predictions. 

The cumulative number of haloes above a given mass Mh is 



N{> Mh) = / dM / 

J M„,ia izi 



dn dV 
dMdV dl 



dz 



(33) 



where dn/dAIdV represents the differential mass function and 
dV/dz the volume element. Since dark energy does not only af- 
fect the growth history but also the geometry, we expect that the 
contribution of volume effects on observable quantities will pro- 
vide more information than merely the differential mass function. 

Effects on the number of observable haloes are shown in 
Fig. [D In the top left panel, we show the ratio of the cumulative 
mass function above a given mass integrated between z — and 
z = 2 between some of the dark energy models studied and the ref- 
erence ACDM model. The top right panel shows the volume effect 
for the corresponding models, i.e. the ratio between the volumes of 
the dark-energy and the ACDM models. In the bottom panels, we 
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Figure 4. The left panels show the time evolution of the linear overdensity 5c{z), the right panels the time evolution for the virial overdensity Ay{z) for 
the different classes of models. In all panels, the ACDM solution (black solid curve) is the reference model, while the black dashed horizontal curve is the 
EdS model that is reached asymptotically by all the models. The upper panels present the quintessence models: the INVl (INV2) model is shown with the 
light-green dashed (dark-green short-dashed) curve, the 2EXP model with the blue dotted curve, the AS model with the cyan dot-dashed curve, the CPL 
(CNR) model with the red dot-dotted (orange dot-short-dashed) curve and finally the SUGRA model with the violet dot-dot-dashed curve. The middle panels 
show the Casimir effect (brown dotted curve) and the (generalised) Chaplygin gas with the (turquoise short-dashed) magenta dashed curve. Finally the lower 
panels report the solution for the models with constant equation-of-state pai'ameter: the dark-green short-dashed curve stands for the non-flat ACDM model, 
the light-green dashed curve for the model with w = —2/3, the blue dotted curve represents the model with w = —4/3, the cyan dot-dashed curve the 
model with w = —1.5, the orange dot-short-dashed curve the model with w = —5/3, the red dot-dotted curve the model with w = —2 and finally the violet 
dot-dot-dashed curve curve shows the model with w = —3. 
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show the contribution to the number counts in spherical shells en- 
closing a volume between z = Q and z = 1 (left panel) and z = 1 
and z = 2 (right panel). Please see the figure caption for details on 
the models considered. 

In the Press & Schechter formalism, the cosmological infor- 
mation, and hence the dark energy contribution, is contained in the 
quantity 5c{z)/ D{z), where D{z) is the growth factor. To com- 
pare the different models, we thus fix the variance for the ACDM 
model (as = 0.8) and scale the variance o f the dark energy models 
according to (see also lAbramo et ahluOOTh 



C"8,DE = 



,{z = Q) 



5c, ACDM (z — 0) 



o-g 



(34) 



Due to the relatively small differences in terms of 5c, the normali- 
sations differ by a few percent at most. 

Even if the differential mass functions differ only slightly, 
we note that differences in the number counts are as large as 
40% — 60%. This is mainly due to volume effects, as shown in 
the top right panel of Fig. [5] Models with a non-phantom equa- 
tion of state always have a smaller volume than the ACDM model 
because the expansion rate is lower, while the opposite holds for 
phantom models (dot-dot-dashed brown curve). Since differences 
in the mass function are expected only in the high-mass tail, we can 
safely assume that the number of small objects is approximately 
the same for all the models. Thus, for low mass haloes, the non- 
phantom dark energy models predict fewer objects, but for objects 
above ~ 10^* Mq/H, the exponential tail of the mass function 
compensates the smaller volume, and we see that a larger number 
of high-mass objects is expected. Of course, for the phantom mod- 
els, the results are reversed. We expect more objects at low mass 
and fewer at high mass, since they do not have time to assemble. 

From an observational point of view, it is also interesting to 
determine in which redshift interval we expect the highest contribu- 
tion. This is shown in the lower panels. We normalise to the ACDM 
counts integrated over the same redshift interval as the dark-energy 
models, thus the sum of the two panels does not reproduce the top 
left panel. It is clear from the bottom right panel that the major con- 
tribution comes from high redshifts, while we do not expect more 
than 10% difference from the volume up to z = 1. Once again, for 
the phantom models, the situation is reversed. 

Despite the fact that differences in number counts are not neg- 
ligible and systematic, we have to recall that they are of the same 
order of magnitude as the uncertainty in the determination of halo 
masses. It will therefore still be difficult to discriminate between 
the models studied. 



6 CONCLUSIONS 

In this work, we have generalised the non-linear equation governing 
the evolution of matter overdensities to cosmological models con- 
taining fluids with an arbitrary equation-of-state parameter. Speci- 
fying w = 0, we recover the well-known equation for structure for- 
mation in matter-dominated universes. By means of the non-linear 
evolution equation, we determine the appropriate initial conditions 
used to solve the linear equation and compute the linear threshold 
for collapse, 5c. We point out that the derivation of Eq. dlSI l is very 
general and can be extended to very broad classes of cosmological 
models, once the appropriate continuity, Euler, Poisson and back- 
ground equations are provided. 

In this work, we considered exclusively non-clustering dark 
energy models, in which the only clustering component is the dark 



matter. Our goal was to study a whole catalogue of dynamical dark- 
energy models, thereby summarising results partly obtained else- 
where, and to clarify discrepant results on early dark-energy mod- 
els. We stuck to the common assumption that the sound speed in 
the dark energy is given by the equation of state. There are more 
general scenarios allowing perturbations also in the dark-energy 
component, and this might lead to importan t differences compa red 
to the case of non-clustering dark energy ( Abramo et aLlbOOTT) . If 
the dark-energy component is not homogeneous, one might won- 
der what happens to the number counts if the equation of state of 
the dark-energy fluid changes within the collapsing sphere. Then, 
the equation of state will depend on the actual overdensity of the 
dark energy, and the model will acquire an additional degree of 
freedom, parametrised by the effective sound speed defined by 
Cs.cff = 5P/5p. Theoretical prediction s for the equation of state 
inside the collapsing sphere are given by A bramo et alj j2008l) . and 
a determination of the effective sound speed in Sunyaev-Zel ' dovic h 
and weak-lensing surveys is carried out in lAbramo et al.l ( 120091) . 
The authors found that a negative pressure perturbation to the dark 
energy fluid may have a substantial effect on the number counts 
as shown in their Fig. 1. A negative effective sound speed c^ can 
lower the value of 5c to ~ 1.5 . . . 1.55, giving a substantial boost 
to structure formation. Volume effects will still be comparable with 
what we found here. It may then be possible to descriminate this 
class of models compared to those studied here. 

Despite the fact that the differential mass functions are very 
similar because of the small differences in 5c, we still expect a sig- 
nificant difference for the total number of objects in a given volume 
above a given mass threshold. For low minimum mass, the volume 
effects dominate and we expect a lower overall number of objects, 
while the mass function dominates over the volume effect at the 
high-mass tail, and we still expect more haloes. Unfortunately, as 
shown in the lower panels of Fig.|5] the major contributions come 
from redshifts above 2 = 1, where fewer objects are expected and 
where observations are more difficult. 

We found that all models studied here show differences at the 
per-cent level compared to the standard ACDM model. We also ar- 
gue that in the framework of general relativity, it may be possible to 
have a more pronounced impact on non-linear structure formation 
if the equation-of-state parameter is discontinuous or if modifica- 
tions to general relativity are involved. This will be the subject of 
future work. 

We thus conclude that, at least for the wide class of models 
studied here and due to the current status of observations, it is 
quite difficult to use number counts to discriminate between dif- 
ferent dark-energy models. A discrimination may still be possible 
to eventually based on geometrical tests where volume effects are 
relevant. 
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Figure 5. Volume effects on lialo number counts. Tlie top left panel shows the cumulative mass function between ^ = and z = 2 normalised with the 
expected value of the ACDM model. The top right panel shows the volume as a function of redshift compared to the volume of a ACDM model. The lower 
panels present the different contributions to the number counts in two redshift bins, between z = and z = 1 and between z = 1 and z = 2in the left and 
right panels, respectively. Different colours refer to different dark energy models. Green short-dashed: INVl; blue dotted: SUGRA; (orange dot-short-dashed) 
cyan dot-dashed: (generalised) Chaplygin gas; red dot-dotted (brown dot-dot-dashed): phantom model with w = —2/3 (w = —3). 
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APPENDIX A: <5c FOR EDE COSMOLOGIES 

Having found that cosmological models with dynamical dark en- 
ergy have very little impact on the spherical-collapse parameters 
Sc and Av, we have to clarify and explain how we could arrive at 
contradictory results in lBartelmann et al.l 1200g). 

We first show the results we obtain with the new approach 
presented in this paper for the linear overdensity 5c in two dif- 
ferent EDE models. We further compare our theoretical results 
to the numerically simulated mass function of iGrossi & Springell 
20091) before we turn t o explain why the earlier calculation of Sc 



Bartelmann et al.l2006r) arrived at significantly different results. 



Figure IaTI shows the time evolution of 5c as a function of red- 
shift (upper panel) as well as the comparison between the theoret- 
ical pred iction for the mass fun ction using the Sheth-Tormen ex- 
pression JSheth & Tormenll 19990 and the numerical mass functions 
obtained via N-body simulations at the redshifts z £ {0, 1, 2, 3} 
(lower panel). 

We notice that EDE cosmologies do not have any strong 
impact on the spherical-collapse parameters and neither on the 
halo mass function. The deviations from ACDM reach at most of 
(1 ... 2)%, in excellent agreement with numerical simulations. 
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Figure Al. Upper panel: Sc for three different early dark energy models. 
The black solid curve represents the ACDM model while the red dashed, 
green short-dashed and blue dotted curves represent the EDEl, EDE2 and 
EDE4 models, respectively. Lower panel: comparison between the theoreti- 
cal multiplicity mass function (given by the differential mass function times 
the mass squared) and the numerical one. Results for four different redshifts 
are shown: 2 = (orange dot-dashed), z = 1 (red dashed), 2 = 2 (green 
short-dashed) and 2 = 3 (blue dotted). Points represent the data of the 
A'^-body simulation. 



In the earlier study by iBartelmann et al.l yOOo), several ap- 
proximations had to be made to render the spherical-collapse 
equations numerically tractable and stable. In contrast to the new 
approach presented here, the main problem there was that the 
spherical-collapse equations become singular at times i — >■ 0, while 
5c must be obtained extrapolating from the limit of the solution 
for i — >■ 0. Yet, the solutions for the radius y{x) of the spherical 
overdensity as a function of the scale factor x as well as for the 
density contrast 5{x), turn out to agree precisely with the solutions 
obtained with the new approach. 

Substantial deviations begin w ith the integration const ant B 
introduced in Eqs. (22) and (23) in lBartelmann et alj ( 120061) . Ac- 
cording to the approximations used there, B « 1 to good accuracy. 
In the new approach, this can be tested using 5 = 3By/5 at early 
times. It turns out that the correct result is B ~ 2.13. This substan- 
tial deviation can be traced back to the integrand (1 + 3w)g{x)y 
of the integral / in Eq. (15). The power-law approximation made 
in Eq. (18) turns out to hold very well for early times, but to fail 
considerably at late times. Thus, the late-time evolution of the in- 
tegrand in Eq. (15) is incorrect, and the deviation from its ex- 
act behaviour starts becoming substantial already at scale factors 
X > 0.1. The power-law approximation of Eq. (19) for the integral 
/ is accordingly incorrect. The assumed power-law behaviour and 
its true shape are shown in Fig. lA2l 
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Figure A2. Exact so lution of the integral / from Eq. (15) in 
iBartelmann et all 120061) compared to the power-law approximation defined 
in Eq. (19) there. The failure of this approximation is the reason for the dis- 
crepant results obtained there compared to the present study. 



Intriguingly, the approximation is correct for conventional 
Friedmann models with arbitrary cosmological constant and also 
for cosmologies with dynamical dark energy with vanishing dark- 
energy density at early times. Its inaccuracy thus remained undis- 
covered in the numerous tests that were carried out. 



